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ABSTRACT 

We introduce the theory of non-linear cosmological perturbations using the correspondence limit of the 
Schrodinger equation. The resulting formalism is equivalent to using the collisionless Boltzman (or Vlasov) equa¬ 
tions which remain valid during the whole evolution, even after shell crossing. Other formulations of perturbation 
theory explicitly break down at shell crossing, e.g. Eulerean perturbation theory, which describes gravitational 
collapse in the fluid limit. This paper lays the groundwork by introducing the new formalism, calculating the 
perturbation theory kernels which form the basis of all subsequent calculations. We also establish the connec¬ 
tion with conventional perturbation theories, by showing that third order tree level results, such as bispectrum, 
skewness, cumulant correlators, three-point function are exactly reproduced in the appropriate expansion of our 
results. We explicitly show that cumulants up to N = 5 predicted by Eulerian perturbation theory for the dark 
matter field 5 are exactly recovered in the corresponding limit. A logarithmic mapping of the field naturally arises 
in the Schrodinger context, which means that tree level perturbation theory translates into (possibly incomplete) 
loop corrections for the conventional perturbation theory. We show that the first loop correction for the variance is 
<t 2 = er 2 +(—1.14+n)cr| j for a field with spectral index n. This yields 1.86 and 0.86 for n = —3, —2 respectively, 
and to be compared with the exact loop order corrections 1.82, and 0.88. Thus our tree-level theory recovers the 
dominant part of first order loop corrections of the conventional theory, while including (partial) loop corrections 
to infinite order in terms of 5. 

Subject headings: cosmic microwave background — cosmology: theory — methods: statistical 


1. INTRODUCTION 


The most successful theories of structure formation assume that small initial fluctuations grew by gravitational amplification in an 
expanding cosmological background. This simple assumption has been remarkably successful in explaining non-linear structures in 
simulations and observations, especially on larger scales. The basic equation governing the statistics of the dark matter field is the 
collisionless Boltzman’s (or Vlasov) equation. These non-linear equations are considered to be intractable, therefore the most usual 
approach is to take their moments. 

Taking only the first few moments of the equation results in Euler’s ideal (pressureless) fluid equations coupled with the Poisson 
equation for gravity. The fl uid approxima tion, however, breaks down at shell crossing. Arbitrarily high moments yield an infinite 
hierarchy of equations (e.g. Peebles, 1980), the BBKGY equations. These relate tim e evolution of fV-th order moments t o N + 1-th 


order moments. The equations can be “closed” only un der certain assumption s (e.g. Davis & Peebles 


1988). Equation of motions in Langrangian space (e.g. Bouchet et al., 1995) are as successful as the Eulerian method. The above 
analytical tools are checked by a set of numerical and semi-numerical methods, such as IV-body simula t ions, and approximations 


1977; Fry, 1984a; Hamilton, 


thereof such as Zeldovich, truncated Zeldovich, frozen flow, adhesion approximation (e.g. Gurbatov et al., 1989). 


One particularly successful method to solve the above equations before shell crossing is perturbation theory (PT); Eulerian (EPT) 
and Langrangian (LPT) perturbation theories are used with si milar success . Cont ributions can be ordered systematically, represented 
and enumerated by the use of Feynman graphs pioneered by poroff et~ak (1986) The first non-trivial order, tree level PT penetrates 


the non-linear regime to surprisingly hig h degree. Cumulants of the dark matter density fie l d have been calculated f or Ga ussian 
initial conditions by several authors (e.g. Peebles, 198()|; Fuszkiewicz et al., 1993|; Bernardeau, 1992; Bernardeau et al., 2002)[]. For 


calculations have been performed for the full three 

-point correlation function and bispectrum (|Ty|. 

1984b 

). The theory have been 

confirmed by simulations (e.g. Colombi et al 


[Szapudi et al., 

19991, 200C; Colombi et al 

, 200C|), and data appear to be in broad 

agreement as well (e.g. Szapudi et alj, 1996; 

Szapudi & Gaztanaga 

1998; Szapudi & Szalay, 

1997; [Szapudi et al. 

2002 

) 


The most straightforward way t o improve tree-level PT is to include the next to leading order, “loop" corrections. Such calculations 
(Scoccimarro & Frieman, 1996b| , jj), albeit fairly complicated, do improve the agreement with simulati ons on small scales. Anothe r 
extension to PT is the spherical infall model, in which one calculates angle averaged Feynman vertices (Fosalba & Gaztanaga, 1998). 


It is a fairly simple alternative to the tedious calculation of loop corrections, but it yields only a fraction of the corrections due to 
neglecting tidal effects. Besides physic al theories, several phen omenological assumptions exist to fit smaller scale behavior (e.g.. 


“extended” and “hyperextended” PT, see Bernardeau et al,, 2002, for details). 


In this paper we propose an entirely different appr oach based on the c orresp ondence limit of the Schrodinger equations. These are 
equivalent to the full Boltzman-Vlasov description (Widrow & Kaiset, 1993), but involve a complex scalar field depending on 3+1 
variables, similarly to the Euler equations. As shown below, perturbation theory of the Shrodiger theory (SPT) is not significantly 
more complicated then EPT. The next section presents the underlying theory, section 3 outlines the connection with conventional PT, 

1 The last paper contains comprehensive references on the subject. 
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section 4 deals with the cumulants in more detail, finally the last section discusses the results and provides an outlook for research 
possible using this formalism. 


2. PERTURBATION THEORY WITH THE SCHRODINGER EQUATION 

The Schrodinger Equation in the correspondence limit is a viable alternative to the collisionless Boltzman’s (Vlasov) equations. 
In the expanding universe (with = 1 for simplicity) 

ihip +—Hip + Hip = 0, (1) 

where t/> is the complex scalar field representing dark matter density with p -\ 'tp | 2 , H is the Hubble’s constant and H is the 
Hamiltonian of the gravitational field. This is a nonlinear equation, since H depends on the density field. If new variables are 
introduced as 


/ \ ~ 3 / 2 

= (2) 

the above equation reduces to equations for two real scalar fields 

^ =-^[V 2 H + 2VHVi?] 

B = ^[V 2 A+ | VA | 2 ] - ^ | VB | 2 -mV 

V 2 V = 4:TrGpa 2 (a 2A — 1), (3) 


where the last equation is the Poisson equation coupled to the Schrodinger equation. The structure of these equations is similar to 
the Euler’s fluid equations in terms of the density contrast and velocity potential, despite the fact that we have not taken moments 
of the underlying full equations. The main difference, and extra complication arises from the exponential in the Poisson equation. 
Perturbation theory of these equations will be entirely analogous to that of the Euler equations, and it can be presented in Fourier 
space in the simplest manner. 

In what follows we will work in the correspondence limit, i.e. h —> 0, we neglect “wavy” features of the equations. Fourier 
transforming the equations yields 


A k — + 2 (kAk)(kBk)\ 

Bk = kB k )(kB h ) - % Ea> r m A k> (4) 

where the multiple and power of transforms are understood as convolutions, m = 1 is assumed for simplicity, the Poisson equation 
was substituted for — V, and the exponential expanded. 

These equations can be rendered homogeneous in a and H using the following Ansatze 


A fc =T,A[ N) a N 

B k = -HY J B[ N) a N + 2 . 


(5) 


Perturbations can be ordered according to powers of the growth factor. We can introduce the usual kernels with the following 
definition 

4 ^ = Jd 3 k 1 ...d 3 k N 6(k = k 1 + ... + k N )F( N \k 1 ,...,k n )A^ ...44 
B l N) = wS ■ • ■ d 3 k N S(k = h+... + k N )G^ N \k u ..., k n )A { y ... 4 . ( 6 ) 

Substituting to the equations leads to the following recursions 

NFW =(?<") +2J2 s a(qi,q2)F^G^ N - s \ 

(N + ^)G (JV) = Es/%r,92 )G( s )G( w - s ) + 3 Em, 8i , 82 ,. = E sJFMfM (7) 

where mode coupling functions are similar to the Eulerian case 


a(<?r,92) 

P(qi,q2) 


_ ( 91 - 92 ) 
^2 ’ 

_ fc 2 (<?l-92) 


( 8 ) 


and the exponential is expanded; < 71 , q%, and k correspond to the sum of S, N — S, and all the wave vectors. To solve the recursion 
at any order, one has to separate F N (M = 1) in the expansion of the exponential and subtract the two equations from each other. 
This procedure leads to terms up to N — 1 on the right hand side of the equations. Here we give explicitly N = 2 case 

3 10 2 

F^(h,k 2 ) = - + —a(k u k 2 ) + ~(3(ki, fc 2 ). (9) 


Higher order kernel functions can be obtained trivially from the recursion relation. These kernels can be used to calculate quantities 
to the accuracy of tree level perturbation theory. In what follows we will show the connection with Eulerian perturbation theory via 
explicitly calculating tree level quantities at third order.. 
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3. CONNECTION WITH EULERIAN PERTURBATION THEORY 

Since S = e 2A — 1, 5\ and S 2 , the perturbative corrections to 6 growing with the first and second power of the growth factor are 


6™ =2 A&, 

S^ = 2(A^+A^f 


( 10 ) 


Thus S 3 to tree level perturbation theory is 

3(<5«V 2 )) 


= 3a 4 8A( 1 r[AW + A( 1 r] 
= 3a 4 8 / d 3 k i... d 3 k 4 ( 1 + F 2 (k 2 , k 3 ))A[\ ] 
= %a 4 fd 3 k 1 ...d 3 k 4 (l + F 2 {k 2 ,k 3 ))s£ . 


-A"e 


lx ki 

E fc *. 


(ii) 


When taking the ensemble average (assuming Gaussian fields), k 2 and k 3 cannot be paired, yielding a combinatorial factor of 2. 
From this S 3 = 3(1 + F 2 ), where () here means angle averaging. Since (a) = 0, and (j3) = 2/3, S 3 = 34/7 as expected from PT. 

In fact it is easy to show 1 + F (2> = 2F <2) (after the necessary symmetrization of the kernels), where f 4 ' 2> is the second order 
kernel (naked vertex) in EPT. Thus calculations entirely analogous to the above show that tree-level bispectrum, 3-pt function, and 
the skewness after smoothing, cumulant correlators, etc. are all exactly matching the tree-level EPT results. 


4. CUMULANTS 


To derive cumulants one has to consider the angle averaged recursion relations. If the angle averaged kernels are defined as 
u/v = N\(F n ) and /zjv = N\{Gn)> the first equation gives a very simple relation: 


Ni/ n = Fn- 

From this the second equation gives a simple recursion for which reads 


v N = 


(21V + 1)N — 3 


A 1 2 / N\ 

£, j s(N - s)v s vNs + 31V! £ ■ 

. a=l d ' s ' M—2 


->M -2 


M\ 


£ S(N = £ ■ ■ ■ — 

^ ^ Si! %! 


( 12 ) 


(13) 


With the initial condition = /xi = 1, one can solve to arbitrary order in a trivial fashion. The first two naked vertices are 

,, — 2 6 
21 

_ 568 
^3 189 

,. _ 473744 
~ 43659 ' 

From these the tree level cumulants of the A field can be calculated as 


(14) 


= 3 ^ 2 = f 

19,,2 _ 40240 

1323 


Si = 4 u 3 

Si = 5124 + 6 OU 3 U 2 + 60^1 = ■ (15) 

Projecting the cumulants of the A field using the formalism of biasing recovers cumulants of 5.. In this context the primary field is 
A and <5 = e 2A — 1 = E bk A k /k\ is a non-line a rly bi ased field, where the usual bias coefficients read 6 = 2 and cn = bN/b = 2 N ~ 4 . 
According to the formula of | 17 & Gaztanaga ( 1993 ), 


S 3 

s 4 

S 5 = b~ 3 (S, 


= 6- 1 (5 3 +3c 2 ) = f, 


= b~ 2 (Si + 


12c 2 Si + 4c 3 + 12c|) = 


60712 


20 c 2 Si + 15c 2 5 : f + (30 c 3 + 120c 2 )^ 


1323 

5c 4 + 60c 3 c 2 ) = 


200575880 

305613 ’ 


(16) 


i.e. we recover the results of tree level PT exactly. 

Because of the non-linear transformation, tree level calculations of the A field translate into (possibly incomplete) loop corrections 
fo r the 6 field. To demonst rate this we calculate the loop corrections arising from the tree-level SPT for the variance a. According 
to Fry & Gaztanaga ( 1993 ) the variance a 2 = 4(A 2 ) + 4(A 2 ) 2 (2Si + 6 ) + .... Initially, there is no skewness, thus for the initial 
conditions, the linear variance ctl follows the same equation with S 3 = 0. Expanding and collecting terms to second order in {A 2 ) 
yields 


2 _ 2 , 

a — a L + 


q A 

D 3 


’L- 


(17) 


Numerically the coefficient is Si/2 = 13/7 ~ 1.857 is within 2% of the exact loop correction 1.82. This means that the dominant 
part of the loop correction arises from tree level perturbations of the logarithm of the field. The same expansion can be generalized to 
any or der, although tedious, fairly simple. It will be present ed elsewhere. Her e we h ave used the general theory of |-ry & Gaztanaga 
(1993), however, exponential bias was explicitly treated by |Grinstein & Wise ( 1986| ); their results will probably be useful for future 
calculations of this sort. 

Note that the above calculations are for the unsmoothed field, or n = —3. Detailed calculation for other spectral indices will be 
shown elsewhere; as a preview, simple considerations suggest Si = 26/7 — 2 (n + 3), and consequently a 2 = a 2 + (—1.14 + n)c 
for n < —1.14 in excellent agreement with PT. 
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5. SUMMARY AND OUTLOOK 


We have introduced SPT based on on the Schrodinger equations. In the correspondence limit, our description is equivalent with 
the full Boltzman (Vlasov) equations. Other versions of PT, in particular EPT, explicitly break down at shell crossing, therefore our 
ultimate aim is to penetrate the non-linear regime deeper. 

A unique feature of SPT that it is naturally uses the logarithm of the dark matter field, which remains close to Gaussian throughout 
the non-linear evolution. This suggests that the SPT expansion should converge faster. Related techniques such as Edgeworth 
expansion, etc. are expected to be more accurate. Indeed, we have shown that, at le ast for the variance, th is is the case. This feature 
of our theory sheds new light on the validity of the lognormal prescription for S (e.g. Coles & Jones, 1991) and suggests using A , i.e. 
the logarithm of the underlying dark matter field to construct statistics, such as correlation functi on, skewness, ku rtosis, bispectrum, 
(in A-space). The first encouraging steps towards analyzing data in log-space have been done by Colombi (1994). 

Our principal aim here is to present the basic theory, and the calculation of the recursions for the tree level perturbation theory 
kernels. These constitute the ground work for a spectrum of future research. We have also given the recursion for the angle averaged 
kernels and elaborated the results up to N = 5 explicitly. Higher orders are calculated trivially from the results shown. 

We have explored the connection of our theory with EPT. We have recovered the first non-trivial tree level PT kernel in the 
appropriate limit. All third order tree level PT results (e.g., bispectrum, three-point function) are exactly reproduced by the SPT 
expansion. In addition, we have applied the exponential bias formalism to recover exactly the tree level PT cumulants up to TV = 5. 
The theoretical exercise of proving the equivalence to arbitrary order is left for subsequent research. 

Our tree level results correspond to (possibly incomplete) infinite order loop corrections. To demonstrate that, we have obtained 
the first non-trivial loop correction for the variance by simply expanding our tree level results to second order. The agreement of this 
simple expansion with the complex EPT loop calculations is remarkable. A large, perhaps dominant, fraction of the loop corrections 
arises from the non-linear projection of tree level results from ,4-space to (5-space. 

The following extensions and generalizations will be presented in subsequent publications: detailed exploration of the higher 
order cumulants, cumulant correlators, A’-point correlation functions and N — 1-spectra, and the probability distribution function; 
other statistics, such as genus and void probability; smoothing; non-Gaussian initial conditions; application to angular correlations 
and lensing; general cosmological background; approximate loop and non-perturbative corrections for higher order quantities using 
our non-linear projection; exact loop corrections in ,4-space; comparison of tree-level statistics of the ,4-field, such as bispectrum, 
three-point correlation function, and cumulants (variance, skewness, and kurtosis), with simulations. Other interesting applications 
of the theory are possible, such as studying the wavy term which was neglected so far (and thus modeling a dark matter with large 
Compton wavelength). A fast approximation scheme similar to 2LPT will be obtained and implemented from our formalism, possibly 
including wavy terms. 

Note that we have not mentioned the interpretation of the B held (G N and //,y). Initially, this corresponds to the usual velocity 
potential; the interpretation after shell crossing is unclear. Nevertheless, our calculations yield G N and /jy as well; the connection 
with measurements is left for future work. 

The method we have put forward follows almost exactly the prescription of EPT for a slightly different set of equations. These is 
not that only way to deal with these equations; other possibilities exist, such as using the interaction picture analogously to quantum 
mechanics. As a first step towards this direction, the authors have shown that the Zeldovich approximation is recovered in the 0-th 
order approximation. This in turn (along with the logarithmic mapping) hints a connection with Langrangian perturbation theory 
(LPT), which will be explored later. 

IS was supported by NASA through AISR grants NAG5-10750, NAG5-11996, and ATP grant NASA NAG5-12101 as well as by 
NSF grants AST02-06243 and PHY99-07949; the latter while enjoying the hospitality of KITP. IS thanks Csaba Balazs for helpful 
discussions. 
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